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Abstract 

o> : 

<N . 

We study the dual Higgs mechanism induced by monopole condensation 

> : 

00 . based on the dual gauge formalism in the maximally abelian (MA) gauge in 

m ■ 

I the lattice QCD. To examine "monopole condensation" in QCD, we study the 

\ monopole part or the monopole-current system appearing in the MA gauge by 

' extracting the dual gluon field B u . First, we investigate the inter- monopole 

potential using the dual Wilson loop in the lattice QCD simulation. In the 
Jjn \ monopole part in the MA gauge, the inter-monopole potential is found to 



X 



be flat, and can be fitted as the Yukawa potential in the infrared region. 



From more detailed analysis of the inter-monopole potential considering the 
monopole size, we estimate the effective dual-gluon mass mg ~ 0.5GeV and 
the effective monopole size Rm — 0.2fm. Second, we study the dual gluon 
propagator G®(x — y) = (B fJi (x)B u (y)}MA in the MA gauge, and find that 
Gj~L behaves as the massive vector-boson propagator with tub — 0.4 GeV in 
the infrared region. The effective-mass acquirement of the dual gluon field 
at the long distance can be regarded as the lattice QCD evidence of "infrared 
monopole condensation" in the MA gauge. 
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I. INTRODUCTION 



Quantum chromo-dynamics (QCD) is the fundamental theory of the strong interaction 
and exhibits the interesting phenomena such as color confinement and dynamical chiral- 
symmetry breaking (D%SB) The QCD Lagrangian has the SU(iV c ) local symmetry 

and is described by the quark field q and the gluon field as 

Cqcd = ~tx(G^ u G^) + q(ip - m q )q, (1) 

where denotes the SU(iV c ) field strength G^ u = -^[D^, D v ] with the covariant-derivative 
operator = + ieA^ Jl|. In QCD, the asymptotic freedom is one of the most important 
features in QCD and the QCD gauge- coupling constant e becomes small in the ultra- 

violet region. On the other hand, in the low-energy region, the gauge-coupling constant e 
becomes large, and there arise the nonperturbative-QCD (NP-QCD) phenomena like color 
confinement and D^SB corresponding to the strong-coupling nature. These NP-QCD phe- 
nomena are extremely difficult to understand in the analytical manner from QCD and have 
been studied by using the effective models |J or the lattice QCD simulation 0. 

In the lattice formalism, space-time coordinates are discretized with the lattice spacing 
a, and the theory is described by the link variable U^i^s) = e iae G SU(iV c ). For instance, 
the standard lattice action is given as 

S = P £ [1 - ^{UM + Ul(s)}}, (2) 

where U^ v (s) is the plaquette variable defined as U^ u {s) = U^{s)U v {s + fj^UUs + u)Ul(s). 
The lattice QCD Monte Carlo simulation is based on the numerical calculation of the QCD 
partition functional, and it is one of the most reliable methods directly from QCD. In the 
lattice QCD, nonperturbative quantities like the quark confinement potential, the chiral 
condensate (qq) and low-lying hadron masses are well reproduced 0. However, one cannot 
understand the confinement mechanism only by looking on the result obtained by the lattice 
QCD. 
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Recently, the lattice QCD simulation has shed light on the confinement mechanism in 
terms of the dual-superconductor picture, which was first proposed by Nambu,'t Hooft and 
Mandelstam in the middle of 1970's |8|-|T0||. In this scenario, quark confinement can be un- 
derstood with the dual version of the superconductivity. In the ordinary superconductor, 
the Meissner effect occurs by condensation of the Cooper-pair with the electric charge [JTTJ . 



Consider existence of the magnetic charges with the opposite sign immersed in the super- 
conductor, then the magnetic flux is squeezed like a tube between the magnetic charges, 
and the magnetic potential between them becomes linear as the result of the Meissner effect 



TT| . On the other hand, the confinement force between the color-electric charge is charac- 
terized by the universal physical quantity of the string tension a = 0.89 ~ 1 GeV/fm, which 
is evaluated from the Regge trajectory of hadrons and the lattice QCD simulations. This 
string tension is brought by one-dimensional squeezing of the color-electric flux in the QCD 



vacuum 



j25| . In the dual-superconductor scenario ||, the QCD vacuum is assumed as the 
dual version of the superconductor, and the dual Meissner effect brings the one-dimensional 
flux tube [JEJ] between the quark and the anti-quark, which leads to the linear confinement 
potential @-§,fnj. 

The dual Higgs mechanism, however, requires "color-magnetic monopole condensation" 
as the dual version of electric condensation in the superconductor, although QCD dose not 
include the color-magnetic monopole as the elementary degrees of freedom. On the appear- 
ance of magnetic monopoles from QCD, 't Hooft showed that QCD is reduced to an abelian 
gauge theory with magnetic monopoles by taking the abelian gauge, which fixes the partial 
gauge symmetry SU(A^ c -)/U(l) 7Vc ~ 1 through the diagonalization of a gauge- dependent vari- 
able [p~4| . Here, the color-magnetic monopole appears as the topological object corresponding 



to the nontrivial homotopy group 7r 2 (S\J (N c ) /\J (l) Nc ~ 1 )=Z^ c ~ 1 . 

Recent lattice QCD studies show abelian dominance ||15|| , which means that NP-QCD 
phenomena like confinement |18|,|19[ and D^BS |20| , |2T| are almost described only by the di- 
agonal gluon component in the maximally abelian (MA) gauge The MA gauge is a 



sort of the abelian gauge where off-diagonal gluon components are minimized by the gauge 



transformation. In the MA gauge, physical information of the gauge configuration is max- 
imally concentrated into the diagonal gluon component. As for the NP-QCD phenomena, 
QCD can be well approximated by the abelian projected QCD (AP-QCD), which is the 
abelian gauge theory extracted from QCD in the MA gauge by removing the off-diagonal 
gluon 0,0. 

AP-QCD includes not only electric currents but also magnetic currents and can 
be decomposed into the photon part and the monopole part corresponding to the separation 
of jn and respectively m other words, AP-QCD consists of the two parts. The 

photon part of AP-QCD (photon-projected QCD) is the abelian gauge theory only with 
electric currents like the ordinary QED. On the other hand, the monopole part of AP-QCD 
(monopole-projected QCD) only includes monopole currents k^. The lattice QCD studies 
show monopole dominance, which means that NP-QCD phenomena are almost described 
only by the monopole projected QCD [p0| , p3| -|27| . In particular, the lattice QCD simulation 
numerically shows that the monopole current is responsible for the electric confinement 
and does not contribute to it in the MA gauge |2^-F25|. Since we are interested in the 
infrared physics such as confinement and its mechanism in the QCD vacuum, we discard 
the off-diagonal gluon and the photon part in the MA gauge as the irrelevant ingredients, 
because they do not contribute to the quark confinement. Then, we concentrate ourselves 
on the monopole part in the MA gauge for the study of the confinement mechanism. 

Then, the remaining problem is whether "monopole condensation" occurs in the QCD 
vacuum or not. In this paper, we perform the SU(2) lattice QCD simulation in the MA gauge, 
and extract the monopole current k^ as the relevant degrees of freedom for confinement. 
Then, we numerically derive the dual gluon field from k^ [p8| ,29[], and investigate the 
effective-mass acquirement of the dual gluon in the QCD vacuum by examining the 
inter-monopole potential and the dual gluon propagator. 
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II. ABELIAN PROJECTION AND MONOPOLE PROJECTION IN THE 

LATTICE QCD FORMALISM 



A. Abelian-Projected QCD in the MA gauge 

Recent studies with the lattice QCD Monte Carlo simulation have revealed the abelian 
dominance and the monopole dominance in the maximally abelian (MA) gauge for the non- 
perturbative QCD (NP-QCD) phenomena such as confinement, dynamical chiral-symmetry 
breaking and instantons fT8|-p7f. In the SU(2) lattice formalism, the MA gauge fixing is 
achieved by maximizing 



R = ltrJ2[T 3 U»( S yul(s)] = E [l - 2 ({^(*)} 2 + {Ul(s)} 2 )} (3) 
by the SU(2) gauge transformation, 

U^Uf A = V(s)U,(s)V\s + ^) } (4) 

where V(s) and V(s + /i) are the gauge functions located at the starting and end points of 
the link variable U^(s). In this gauge, the absolute value of off-diagonal components U^s) 
and U^(s) are forced to be small as possible using the gauge degrees of freedom. 

In accordance with the Cartan decomposition, the SU(2) link variable U^s) is factorized 



as 



U^ A (s) = M,(s)u,(s) , (5) 
M^exp^^ + ^T 2 )], u M ( S )=exp[z(^r 3 )], (6) 

where w M G U(l) 3 and M M G SU(2)/U(1) 3 correspond to the diagonal part and the off- 
diagonal part of the gluon field, respectively. In the continuum limit, the angle variable Q a 
goes to the gluon field A a ^ as Q a — > \eaA a ^. The off-diagonal factor M M (s) is rewritten as 



MM = e^ Tl+e ^ 



( cos#„ — sin#„e 
ysin^e iXM cos 9^ j 



(7) 



with 







^ = mod^^^ + {6l)\ X , = tan- 1 ^. (8) 

Here, the parameter ranges are usually taken as 

^e[0,|], [0,2tt), ^e[0,27r), (9) 

which provides the one-to-one correspondence to the SU(2) group element. Under the abelian 
gauge transformation with v (s) e U(l)3, M M (s) and u^s) are transformed as 

M,(s) - M;( S ) = v(s)M,(s)v\s), (10) 

-> U*(s) = V{8)u li (s)v^[s + /i), (11) 

to keep the form of Eq.(§) for e SU(2)/U(1) 3 and e U(l) 3 . Then, M M (s) behaves as 
the charged matter field and the abelian link- variable 



( e ifl» o ^ 
i e"^ i 



(12) 



behaves as a abelian gauge field with respect to the residual abelian gauge symmetry. 

As a remarkable feature of the MA gauge, the abelian dominance holds for the NP- 
QCD phenomena such as quark confinement and chiral- symmetry breaking [pj|-|2~T|| . Here, 
we call abelian dominance for an operator O, when the expectation value (0[I7 M ]) is almost 
equal to the expectation value (0[uJ}ma, where off-diagonal gluons are dropped off in the 
MA gauge. For instance, the abelian string tension OAbei = (^(^Dma in the MA gauge 
is almost equal to o"su(2) = ( cr (^)) as ^Abci — 0.92cr SU (2) for (3 ~ 2.5 in the lattice QCD 
53J2S[. Thus, NP-QCD phenomena are almost reproduced only by the abelian link variable 



Up, and off-diagonal gluon components 9 l , 9 2 do not contribute to NP-QCD in the MA 
gauge. Hence, as long as the infrared physics is concerned, QCD in the MA gauge can be 
approximated by the abelian projected QCD (AP-QCD), where the SU(2) gluon field 8^T a 
is replaced by the abelian gluon field O^t 3 . In other words, AP-QCD is the abelian gauge 
theory keeping essence of NP-QCD, and is extracted from QCD in the MA gauge. Hereafter, 
we pay attention to the AP-QCD described by 8 3 in the MA gauge. 



B. Decomposition of AP-QCD into Monopole Part and Photon Part 



The abelian-projected QCD (AP-QCD) includes not only the color-electric current 
but also the color- magnetic monopole current [14]. In the Maxwell equations with j M and 
kfj,, the field strength 9^ el satisfies as 

^C* = 3u (13) 
d, *C el = k » ( 14 ) 

with *0^ el = ^e^ upa 9p^ el . In the presence of both j M and k^, the field strength 9^ el cannot 
be described by the simple two-form dffi — d u 9^ with the regular one-form 9^ [0O[| . Using 
the diagonal gluon component 0^, the abelian field strength 9^ el is defined by 

9^ cl = mod^M - d v 9l) G [-7T, tt), (15) 

which is U(l)3-gauge invariant. In the continuum limit a — ► 0, 9^ eX goes to the abelian 
field strength F^ u in the physical unit as 9^ cl — > ^ea 2 ^^, and therefore the range of 9^ cl 
is required to include 0. Then, the two- form 8^91 — d u 9^ can be decomposed as 

dX-d„9l = 9^ el + 27rn^, (16) 

where the former part denotes the field strength and the latter part 27rn^(s) G 27rZ cor- 
responds to the Dirac string on the lattice |22[ . Thus, in the ordinary description |30| . 



the system includes the singularity as the Dirac string 2irn^ u , which makes the analysis 
complicated. 

To clarify the roles of and k^ to the nonperturbative quantities of QCD, we consider 
the decomposition of AP-QCD into the photon part and the monopole part, corresponding 
to the separation of and k^. We call this separation into the photon and monopole parts 
as the "photon projection" and the "monopole projection", respectively p6| , |3T| . 

In the lattice formalism, the photon part 9^ h (s) and the monopole part 9jf°(s) are ob- 
tained from 9f^ el (s) and 27m Atl ,(s), respectively, 



(17) 
(18) 



using the inverse d'Alembertian □ 1 on the lattice [20,32|- Here, □ 1 is the non-local 



operator p^j34 



1 1 



47T 2 (s - S') 



l\2 



(19) 



which satisfies 



U s ( s \U- l \s') = 5\s - s') = (sis'). 



The diagonal gluon component 9^(s) is found to be decomposed as 



(20) 



iMo/ 



(21) 



in the Landau gauge, d^9^(s) = 0. The field strengths, 9f^ in the photon part and 9^° in 
the monopole part, are given as 



C = mod 27r (aX h - dj 



Ph 



9 Ph ) 



= mod 2 ,(^ Mo -^ 



Mo\ 



G [-7T,7r) 
G [ 7T, 7T*) 



(22) 
(23) 



on the lattice, respectively. 

In the actual lattice QCD simulation, the monopole current fc M and the electric current 
are slightly modified through the monopole and the photon projections, respectively, due 
to the numerical error on the lattice. However, these differences are negligibly small in the 



actual lattice QCD simulation. In fact, k^° ~ k^ and ~ hold in the monopole part, 
and j Ph ~ jfj, and k^ h ~ hold in the photon part within 1% error ||26|| . Here, we have kept 
the labels as "Mo" and u Ph" for the electric current and the monopole current, and we 
have used (k^°,j^f°) and (A£ , ) for these currents in the monopole part and the photon 
part, respectively. 

As a remarkable fact, lattice QCD simulations show that nonperturbative quantities such 
as the string tension, the chiral condensate and instantons are almost reproduced only by 
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the monopole part in the MA gauge, which is called as monopole dominance p3|-|26|. On the 
other hand, the photon part dose not contribute these nonperturbative quantities in QCD. 

Since we are interested in the NP-QCD phenomena, it is convenient and transparent to 
extract the relevant degrees of freedom for NP-QCD by removing irrelevant degrees of freedom 
like the off- diagonal gluons (6 1 , 6 2 ) and the electric current j M . Therefore, we concentrate 
ourselves to the monopole part, which keeps the essence of NP-QCD such as confinement, 
and we examine whether monopole condensation occurs or not. In other words, our aim is 
to investigate the feature of the monopole- current system appearing in the MA gauge. 



10 



III. DUAL GAUGE FORMALISM 



In the MA gauge, the monopole part carries essence of the nonperturbative QCD like the 
electric confinement |20,23-p?jj. Since the monopole part only includes the color-magnetic 



current k^, i.e. j^—O, the Maxwell equation in the monopole part becomes 

M 0J» = (24) 
d/C = K, (25) 

where 9^° denotes the field strength in the monopole part. This system resembles the dual 
version of QED with ^ and k^ = 0, and hence it is useful to introduce the dual gluon 
field B^(s) instead of 9^°(s) in the monopole part for the study of the dual Higgs mechanism 
in QCD H|. 

The dual gluon field is defined so as to satisfy the definition of the abelian gauge 
field, 

d„B„-d„B lx =*9f°, (26) 

which is the dual version of the ordinary relation, F^ u = d^A u — d u A^ p9 |. The interchange 
between and B^ corresponds to the electro-magnetic duality transformation, H <-> E. 
Therefore, owing to the absence of j^, the dual gauge field B^ can be introduced without the 
singularity like the Dirac string. In other words, the absence of j u is automatically derived 
as the dual Bianchi identity, 

> = %, = 9/(9A% = 0. (27) 

Let us consider the derivation of the dual gauge field B^ from the monopole current k^. 
Taking the dual Landau gauge d^B^ = 0, the relation *9^° = d 2 B u — d^d^B^) = k u 
becomes a simple form OB^ = k^. Therefore, the dual gluon field B^(x) is obtained as 

B v {x) = (n-%) (x) = I d 4 y(x\a-'\y)k u (y)- = / ^T^Z^a ( 28 ) 
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by using the inverse d'Alembertian □ 1 [p2,p2|. Here, the dual gauge formalism provides 



the natural description of the monopole part in terms of the monopole current k^ and the 
dual gluon B^. Moreover, the dual gauge formalism is useful to examine the dual Higgs 
mechanism in QCD [ f28| , [2"9|1 . 

In the dual superconductor picture in QCD, and correspond to the Cooper-pair 
and the photon in the superconductor, respectively. The Cooper-pair and the photon are 
essential degrees of freedom which bring the superconductivity. In the superconductor, the 
photon field gets the effective mass as the result of Cooper-pair condensation, and this 
leads to the Meissner effect. Accordingly, the potential between the static electric charges 
becomes the Yukawa potential Vy(r) oc in the ideal superconductor obeying the London 
equation. Similarly, the dual gluon B^ is expected to be massive in the the monopole- 
condensed system, and the mass acquirement of leads to the dual Meissner effect. In 
other words, the acquirement of dual gluon mass m# reflects monopole condensation, and 
brings electric confinement. Hence, we can investigate the dual Higgs mechanism in QCD by 
evaluating the dual gluon mass m#, which is evaluated from the inter- monopole potential. 

To estimate the inter-monopole potential, we propose the dual Wilson loop Wd |2"E| , |2"9"|] . 
The dual Wilson loop Wo is defined by the line-integral of the dual gluon field along a 
closed loop C, 

W D (C) = i?e[exp (i I B^dx^l (29) 
J c 

which is the dual version of the abelian Wilson loop Wp^ e \{C) = i?e[exp (i § c 8^dx^)}. Here, 
the dual Wilson loop Wd(R x T) describes the interaction between the monopole-pair with 
the test magnetic charges | and — |. The inter-monopole potential^ is obtained from the 
dual Wilson loop as 

V M (R) = - Km ±-\n{W D {R x T)> (30) 



1 Rigorously, Vm(?") corresponds to the inter-monopole potential in the monopole part (the 
monopole-current system) in the MA gauge. 
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in a similar manner to the extraction of the inter-quark potential from the Wilson loop @-[7j . 

To summarize, starting from the monopole current in the MA gauge, we have intro- 
duced the dual gluon field and the dual Wilson loop for the investigation of the dual 
Higgs mechanism in QCD. 
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IV. INTER-MONOPOLE POTENTIAL AND DUAL GLUON EFFECTIVE MASS 



In this section, we show the numerical result of the lattice QCD simulation. For the 
study of the dual Higgs mechanism in QCD, we calculate the dual Wilson loop Wd{R,T) 
and investigate the inter-monopole potential Vm{t) in the MA gauge using the SU(2) lat- 
tice with 20 4 and (3 = 2.2 ~ 2.3. All measurements are performed at every 100 sweeps 
after a thermalization of 5000 sweeps using the heat-bath algorithm. We prepare 100 sam- 
ples of gauge configurations. These simulations have been performed using NEC SX-4 at 
Computation Center of Osaka University. 

In the dual Higgs mechanism, it is essential that the dual gluon field acquires the 
effective mass. For investigating the effective-mass acquirement of the dual gluon which 
is brought by the monopole condensation, let us extract the inter-monopole potential Vm{t) 
from the dual Wilson loop (Wd(R,T))ma obtained by the lattice QCD. The dual Wilson 
loop (Wd(R, T))ma seems to obey the perimeter law rather than the area law for large loops 
as shown in Fig.[|. Since the dual Wilson loop (Wd{R,T))ma satisfies the perimeter law as 

\n(W D (R x T)) MA a -2a(R + T) ■ a (31) 

for large R and T, the inter-monopole potential becomes constant 2a in the infinite limit of 

T, 

2a 

Vm(R) -> lim — R + 2a = 2a. (32) 

T^oo T 

However, in the actual lattice calculation, we have to take a finite length of T, and hence, the 
linear part (2a/T)R remains as a lattice artifact. Therefore, it is necessary to subtract this 
lattice artifact {2a/T)R for evaluating of the inter-monopole potential Vm{t) from Wd{R, T) 
in the lattice QCD simulation. The parameter a can be estimated from the slope of the 
dual Wilson loop \ti(Wd(R X T))ma for large R and T for each lattices. Here, we obtain 
a ~ 0.14GeV from the dual Wilson loop in Fig.|l| for R, T > 3 in the lattice unit. After 
the subtraction of the lattice artifact (2a/T)R, we consider the shape of the inter-monopole 
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potential Vm{t) in the lattice QCD in the physical unit with r = aR. The inter-monopole 
potential Vjw(r) is short-ranged and flat as shown in Fig.Q. 

Here, we compare the inter-monopole potential Vm{t) with the Yukawa potential 

Mr) = (33) 

Since the Yukawa potential Vy(r) satisfies the relation ln{rVy(r)} = -m^r + const, we 
show the logarithm plot of rl/M ( r ) a s the function of r in Fig.|3|. In the long-distance region, 
ln{rVjv/(^)} seems to decrease linearly with r. From the linear slope of ln{rVM(?")} m the 
region of r > 0.4 fm, the effective mass of the dual gluon is estimated as m# ~ 0.5 GeV. 
As shown in Fig.|], the inter-monopole potential can be fitted by the Yukawa potential Vy(r) 
in the long distance region using the dual gluon mass tub = 0.5GeV and the gauge coupling 
e = 2.5. 

Finally, we consider the possibility of the monopole size effect, because the monopole is 
expected to be a soliton like object composed of gluons. In fact, from the recent lattice QCD 
study, QCD monopoles include large off-diagonal gluons component in their internal region 
even in the MA gauge . We introduce the effective size R M of the QCD-monopole, 

and assume the Gaussian-type distribution of the magnetic charge around its center, 



1 /-Ix 



2' 



p ^ RM) ^Tv^r exj) vw l (34) 



Since the monopole part is an abelian system, simple superposition principle on is ap- 
plicable like the Maxwell equation. Therefore, the Yukawa-type potential Vm{t) with the 
effective size R M of the monopole is expected to be 

\r< d ^ ( e / 2 ) 2 f J* ( » W u xexp(-m B |x- Xl +x 2 |) 
V(x]R M ) = / d xi / d s 2 p(xi; i? M )p(x 2 ; R M ) ; ; , (35) 

AlT J J |X — X1+X2I 

or equivalent ly 

V(r;R M ) = [°° dr x f°° dr 2 e- {r " +r ^ /R M f d6 x f dB 2 sin 6 X sin 6 2 

n Rm Jo Jo Jo Jo 



exp 
x 



mB \ { V ( r — r2 cos ^2) 2 + (j"2 sin 9 2 ) 2 — T\ cos^i} 2 + (ri sin 6*i) 2 



(36) 



{J(r — r 2 cos^) 2 + (7*2 sin 6 2 ) 2 — T\ cos^i} 2 + {v\ sin 6*i) 2 
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where r = |x — y| is the distance between the two monopole centers. We apply this poten- 
tial V(r;R M ) to the inter-monopole potential Vm(V) in Fig^. The Yukawa-type potential 
V(r;R M ) with the effective monopole size R M = 0.21fm seems to fit the inter-monopole 
potential well in the whole region of the distance r. 

Thus, we estimate the dual gluon mass mg ~ 0.5GeV and the effective monopole 
size R M — 0.2fm by evaluating the inter-monopole potential Vm{t) from Wd(R,T) in the 
monopole part in the MA gauge. We find the effective mass acquirement of the dual gluon, 
which is essential for the dual Higgs mechanism in the dual superconductor scenario, and 
suggest that color confinement is caused by the monopole condensation in the QCD vac- 
uum. The monopole size R M would provide the critical scale for the nonperturbative QCD 
in terms of the dual Higgs theory, because off-diagonal gluons contribute to the physics and 
AP-QCD should be modified at the shorter scale than R M , as well as the structure of the 
't Hooft-Polyakov monopole 0. 



16 



V. DUAL GLUON PROPAGATOR AND DUAL GLUON EFFECTIVE MASS 



In the previous section, we estimate the dual gluon mass m# using the inter- monopole 
potential. The more detailed measurement of from the inter-monopole potential seems 
difficult because it requires the twice logarithmic-part procedure from the dual Wilson loop 
25fl. As alternative method, we investigate the dual gluon mass mg through the dual 



gluon propagator in the lattice QCD. 



A. Massive Gauge Field and its Propagator 

In this subsection, we consider the general argument on the analytical relation between 
the massive gauge field and its propagator. To this end, we idealize here the dual gluon field 
as a massive vector boson. 

The lagrangian for the free massive vector field in the Proca formalism is written as 

£proca = ~\{d^B u - d v Brf - l -m%Bl (37) 

in the Euclidean metric [|]]. The Proca lagrangian Eq.(|3"TD leads to the field equation 

d^By - d v B^) - m%B v = 0. (38) 

Taking the divergence of this equation, we find 

mldpBn = (39) 

For the massive case m# ^ 0, B^ is divergenceless as d^B^ = and Eq.(|38|) reduces to the 
Klein-Gordon equation 

(□-m!)B„ = 0. (40) 

The vanishing of d^B^ means that one of the four degrees of freedom of _B M is eliminated 
in a covariant way. In this case, the propagator G® u (k) of the massive vector boson B^ is 
derived as 
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in the momentum representation [[|. The propagator G® u (x) in the coordinate space is 
obtained by performing the Fourier transformation as 

G D (x-y) = (B,(x)B u (y)} = [ -^J k < x ^G D (k) 



d k ;i.t„_„,\ 1 ( r k,,k 



(27r) 4 k 2 + m B \ m B J 



Here, we consider the scalar-type propagator, 



_ e ik-(x-y) ( 3 + 1 » 



(2tt) 4 yP + ml m|y 

* e tt-(«-v) _i_ — _ — §^(x — v) (43) 

(2tt) 4 fc 2 + m| m| K VU 1 ' 

because (x — y) depends only on the four- dimensional distance r E = \x — y\ and it is 
convenient to examine G^ (r E ) for estimating the mass of the field B^. The integration in 
Eq.(|43|) is found to be expressed with the modified Bessel function Ki(z), 

2 /0 °rfA: fc ! Sin2 ! ( r d6 sin 2 Be** cosd 



(27r) 3 Jo fc 2 + m% \Jo 



7rJ jo k z + m B \J-i J (2tt) 2 Jo k z + m B 

1 m„ 

Ki(m B rJ. (44) 



(2vr) 2 r f 

Here, we use the integration formula for the Bessel function J\{z 



J x {z) = J^dtVl^¥e izt , (45) 



and the formula for the modified Bessel function, 



00 



2 



zK x {z) = / dt- r — i J 1 {t). (46) 



t 2 + z 2 



Thus, the scalar-type propagator G^Jr B ) is written as 



Cfe) = ^^KMBr E ) + -V(x). (47) 



Using the asymptotic expansion 



- £ Jm-"»)(w.)- (48) 



the scalar-type propagator Eq.(|4"7|) reduces to 



«?„fe) - l#S^ (49) 



wVe ' 2(2tt)- 2 ,-■ 



at the long-range as r £ > m^ 1 L , since the modified Bessel function satisfies the relation 
^Ki(z) ~ ^J^z~ 3 / 2 e~ z for z > 1 as shown in Fig|| In the Eq. (^9|) , the damping factor 
e -m B r E eX p resses short-range interaction mediated by the massive field in the coordinate 
space. Then, the mass mg of the vector field B^x) is estimated from the slope in the 

3 

logarithmic plot of -^-G^{r E ) as the function of r E , 

In J -^G^(rJ \ - -m B r E + const. (50) 

[ V m B ) 

B. Numerical Results from Lattice QCD Study 

In this subsection, we show the numerical results of the lattice QCD simulation on 
the propagator and the effective mass of the dual gluon B^. For the study of the dual 
Higgs mechanism in QCD, we numerically calculate the dual gluon propagator G^ (r E ) = 
(B li (x)B l j,(x))yLA and try to estimate of the effective mass ttib of the dual gluon field B^ in 
the MA gauge using the SU(2) lattice with 24 4 and (3 = 2.4, 2.45. All measurements are 
performed at every 200 sweeps after a thermalization of 10000 sweeps using the heat-bath 
algorithm. We prepare 60 samples of gauge configurations. These simulations have been 
performed using NEC SX-4 at Computation Center of Osaka University. 

For the effective-mass acquirement of the dual gluon B^ as the result of infrared monopole 
condensation, we evaluate the dual gluon propagator obtained from the lattice QCD in the 
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MA gauge. Since the massive vector-boson propagator is expected to behave as G^Jr E ) ~ 
r~ 3 / 2 exp(— rriBr E ) at the long distance as shown in Eq. (|50|) , we investigate the logarithm 
plot of r E 3 ^ 2 G^ (r E ) as the function of r E for the estimation of the dual gluon mass m^. 

We show in Figs||0 the lattice QCD data for the dual gluon correlations, hi{r E ^ 2 G^{r E )} 
and G^(r E ), as the function of the 4-dimensional Euclidean distance r E . In Fig.[| the dual 
gluon correlation ^{f E ^ 2 G^ (r E )} decreases linearly with r E in the region of r E > 0.8 fm. 
We try to fit a straight line to the lattice QCD data of ^{r E ^ 2 G^(r E )} in this region, and 
estimate the effective mass mg of the dual gluon as m# ~ 0.4 GeV (tub — 0.41 GeV 
for (3 = 2.4 and mg — 0.38 GeV for (3 = 2.45). In Fig.|6](b), based on Eq . (|50"D , we compare 
the dual gluon correlation In r E )) with the massive vector-boson correlation with 

rriB = 0.4 GeV in the long distance as r £ ^ 0.8 fm. We find good agreement between the 
lattice QCD data and the massive boson case denoted by the dotted line, and then the dual 
gluon field is consider to have the effective mass itlb — 0.4 GeV in the infrared region. 

In Fig.0, we show the scalar- type dual gluon propagator G®(r E ) = (B^(x)B^(y))uA in 
the linear scale. The dotted curve denotes the analytical form in Eq.(^) for the massive 
vector-boson propagator with m# = 0.4 GeV. In the long-distance region as r E >0.8 fm, we 
find 

G ^{r E ) ^ const ■ —Ki{m B r E ), (51) 

and the dual gluon B^ seems to propagate as the massive vector field with tub — 0.4 GeV. 
On the other hand, in the short- distance region, G®^(r E ) seems to differ from the massive 
vector propagator. We speculate that such a short-distance deviation on G^(r E ) is brought 
by the monopole-size effect, because the monopole is considered as a soliton-like object 
composed of gluons P5|| . 

In this way, we find the effective mass mg ~ 0.4 GeV of the dual gluon B^ from the 
infrared behavior of the dual gluon propagator G^ u = {B^^B^y^MA in the MA gauge 
using the lattice QCD simulation. The effective-mass acquirement of the dual gluon B^ 
supports the dual Higgs mechanism by monopole condensation in the QCD vacuum. 
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VI. SUMMARY AND CONCLUDING REMARKS 



In order to clarify the dual superconductor picture for the quark confinement mechanism 
in the QCD vacuum, we have studied the effective-mass acquirement of the dual gluon field, 
which is essential for the dual Higgs mechanism. In the maximally abelian (MA) gauge, 
QCD is reduced into an abelian gauge theory including both the color-electric current and 
the color-magnetic monopole current. 

For the investigation of the dual Higgs mechanism by monopole condensation in QCD, 
we have introduced the dual gluon field and have studied its features in the monopole 
part (the monopole-current system) using the MA gauge in the lattice QCD Monte Carlo 
simulation. Owing to the absence of the electric current, the monopole part resembles the 
dual version of QED, and hence this system is naturally described by the dual gluon field 
without meeting the difficulty on the Dirac-string singularity. In the dual gauge formalism, 
the dual Higgs mechanism is characterized by the acquirement of the effective mass of 
the dual gluon field B^. Then, we have investigated the dual gluon mass by examining 
the inter-monopole potential and the dual gluon propagator, respectively. 

As the first attempt to evaluate the dual gluon mass, we have calculated the dual Wilson 
loop (Wd(R x T))ma, and have studied the inter-monopole potential Vm{t) in the monopole 
part in the MA gauge by using the lattice QCD simulation. In the lattice QCD, we have 
found that the dual Wilson loop obeys the perimeter law for large loops. Considering the 
finite-size effect of the dual Wilson loop, we have numerically derived the inter-monopole 
potential Vm(^), and have found that Vm{t) is short-ranged and flat in comparison with the 
linear inter-quark potential. Then, we have compared the inter-monopole potential Vm{t) 
with the Yukawa potential and have estimated the dual gluon mass as mg ~ 0.5GeV, which 
is consistent with the phenomenological parameter fitting in the dual Ginzburg-Landau 
theory fL3| . The generation of the dual gluon mass m# supports the realization of the dual 
Higgs mechanism and monopole condensation in the QCD vacuum. To explain the short- 
range deviation between the inter-monopole potential Vm{t) and the Yukawa potential, we 
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have considered the effective size R M of the monopole, which is considered as a soliton-like 
object composed of gluons. We have found a good agreement of Vm{t) and the Yukawa- 
type potential V(r; R M ) with the monopole effective size R M in the whole region of r. The 
monopole size has been estimated as R M — 0.2fm, which would provide the critical scale for 
the NP-QCD in terms of the dual Higgs theory, because the monopole should be treated as 
a non-local soliton at the shorter scale than R M ■ 

As the second attempt, we have evaluated the effective dual-gluon mass by investigating 
the dual gluon propagator G^(x — y) = (B fJ/ (x)B fJi (x))uA in the MA gauge in the lattice 
QCD. The dual gluon propagator seems to be well fitted by the massive vector-boson prop- 
agator at the long distance. From the behavior of the dual gluon propagator in the infrared 
region, we have estimated the effective dual-gluon mass m# ~ 0.4 GeV. 

Thus, we have shown the effective mass acquirement of the dual gluon field in the 
studies of the inter-monopole potential and the dual gluon propagator in the MA gauge 
using the lattice QCD simulation. This result support the dual Higgs mechanism induced 
by the monopole condensation, which would be responsible for quark confinement, at the 
infrared scale of QCD. The values of the dual gluon mass obtained by the two methods seem 
to coincide within the numerical error. Here, m^ 1 ~ 0.4 fm corresponds to the "penetration 
depth" of the dual superconductor and provides the radius of the color-electric flux tube of 



hadrons ||13|| . Then, the estimation of the dual gluon mass, mg = 0.4 ~ 0.5 GeV, determines 
the most relevant quantity of the dual superconducting theory or the key parameter of the 
dual Ginzburg Landau (DGL) theory |13|]. 
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FIGURES 



FIG. 1. The dual Wilson loop (W D (R,T)) MA as the function of its perimeter L = 2(R + T) 
in the monopole part in the MA gauge in the SU(2) lattice QCD with 20 4 lattice and j3 = 2.3. 
The data for the large loops with R,T > 3 are shown. The perimeter law seems to hold for 
(W d (R,T))ma. 

FIG. 2. The inter-monopole potential Vm{t) as the function of the distance r between the 
monopoles in the monopole part in the MA gauge in the SU(2) lattice QCD with 20 4 lattice and 
P = 2.207,2.257,2.3. 

FIG. 3. The logarithmic plot for tVm(^) as the function of r in the monopole part in the 
MA gauge in the SU(2) lattice QCD with 20 4 and (3 = 2.207,2.257,2.3. In the long-distance 
region, ln{rVM(r)} seems to decrease linearly with r. For comparison, We plot also the solid line 
ln{rVM(r)} = —msr + const with tub = 0.5 GeV. The effective mass of the dual gluon is 
estimated as tub — 0.5 GeV from the linear slope of ln{rVM(r)} in the region of r > 0.4 fm. 

FIG. 4. The inter-monopole potential Vm{t) as the function of the distance r = |x — y| between 
the monopoles in the SU(2) lattice QCD with 20 4 and (3 = 2.207, 2.257, 2.3. The dotted and the 
solid curves denote the simple Yukawa potential Vy(r) and the Yukawa-type potential V(r;R M ) 
with the monopole size effect, respectively. The Yukawa-type potential V{r;Rm) seems to fit the 
inter-monopole potential Vm^t) with the effective monopole size Rm ^ 0.2fm in the whole region 



FIG. 5. (a) The function G(z) = -K\{z) (solid curve) appearing in the massive boson propa- 
gator and its asymptotic form J j^e~ z (dotted curve) as the function of z. (b) The logarithmic 




of r. 
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FIG. 6. (a) The logarithmic plot of ^ 3/,2 C^(r^ ) as the function of r E in the monopole part in the 
MA gauge in the SU(2) lattice QCD with 24 4 and (5 = 2.4,2.45. The correlation In (r* /2 G^(r E )) 
decreases linearly with r E in the long-distance region, (b) The comparison with the analytical 
form in Eq.(|47|) for the massive vector-boson propagator with the mass mg = 0.4 GeV. The 
massive-boson correlation denoted by the dotted line almost reproduces the lattice QCD data for 
the dual gluon correlation In (r E ^ 2 G^(r E )^j at the long distance as r ^ 0.8 fm. The effective mass 
of the dual gluon is estimated as tub — 0.4 GeV from the linear slope of In (j" E ^ 2 G^(r E )j in the 
infrared region. 

FIG. 7. The scalar-type dual gluon propagator G^(r E ) as the function of the 4-dimensional 
Euclidean distance r E in the monopole part in the MA gauge in the SU(2) lattice QCD with 24 4 
and (3 = 2.4, 2.45. The dotted curve denotes the propagator of the massive vector boson with the 
mass tub = 0.4 GeV. The dual gluon seems to propagate as the massive vector field at the long 
distance as r E > 0.8 fm. 
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Figure g 
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